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Advances in quantum technology increasingly 
rely on combining the functions of different sys- 
tems. Strongly interacting systems, such as ions 
[U-tH or superconductors [S], are ideal for process- 
ing; large ensembles for memory Q|; and optical 
photons for communication [5]. However, inter- 
facing these components remains a challenge. For 
example, although cavity QED in the microwave 
domain, using Rydberg atoms [6J or supercon- 
ducting circuits [7j, provides efficient coupling be- 
tween photons and static qubits, microwave pho- 
tons are not ideal for quantum communication 
due to the blackbody background. For this rea- 
son, quantum interfaces that combine different 
functions of a network are desirable. Here we 
demonstrate an interface between the microwave 
and optical domains that allows processing of op- 
tical photons using microwave fields 8|. We store 
optical photons in highly excited collective states 
(Rydberg polaritons) of a cold atomic ensemble 
|9l4ll|]. Subsequently, we apply a microwave field 
that performs a rotation of the collective spin 
of the Rydberg polaritons, leading to coherent 
many-body Rabi oscillations in the amplitude of 
the retrieved photon pulse. We show that the mi- 
crowave field modifies the long-range interactions 
between the stored photons, a key step toward 
the realisation of an all optical analogue of neu- 
tral ato m quantum gates based on dipole blockade 
(l2L Il3l |. In addition, Rydberg polaritons provide 
a powerful platform for studying strongly coupled 
atom light interactions without a cavity, quan- 
tum simulation of resonant energy transfer |l4| 
or spin liquids Il5ll. and quantum metrology us- 
ing Dicke states [16]. 

To convert optical photons into Rydberg polaritons, 
we use electromagnetically induced transparency (EIT) 
pjj ■ A strong control field couples photons in an optical 
signal field to a highly excited Rydberg state, with prin- 
cipal quantum number n [l8|. By reducing the control 
field intensity, the signal pulse is compressed and stored 
in a coherent superposition of the ground and Rydberg 
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states. This collective state is known as a dark-state po- 
lariton [l9| . The dipole-dipole interaction between Ry- 
dberg atoms requires that the polaritons are separated 
by a distance R > R a = (C 6 /?iAeit) 1/6 [IS, known as 
the "blockade radius" for optical excitation, where Cq 
is the van der Waals coefficient that scales as n 11 and 
Aeit is the EIT linewidth. Thus by setting the principal 
quantum number n, and controlling the geometry of the 
ensemble, one can write a desired number of polaritons, 
M. 

We perform quantum state control of the polaritons us- 
ing a microwave field resonant with a close-lying Rydberg 
state. Under the action of this dressing field, each polari- 
ton behaves as a collective two-level system (see Supple- 
mentary Information). The microwave field also modifies 
the dipole-dipole interaction between polaritons [U [22| , 
and introduces a second blockade scale with a character- 
istic size i? M = (Ca/M^) 1 ' 3 , where C3 is the resonant 
dipole-dipole interaction coefficient and f2 M is the Rabi 
frequency of the microwave transition. By varying fi M 
we can tune the ratio R^/Ro- For i? M < R Q the polari- 
tons arc effectively non-interacting and one observes M- 
particle Rabi oscillations that couple to a single optical 
radiation mode. In contrast, for R^ > R a the polaritons 
are strongly interacting and the dynamics are dominated 
by resonant energy exchange [3, HH • Consequently an 
array of polaritons provides an analogue of a spin liquid 
with tunable nearest-neighbour interactions. 

A schematic of the experiment is shown in Fig. 1. Fur- 
ther details are given in the Methods Summary. An op- 
tical dipole trap confines an ellipsoidal atomic cloud con- 
taining up to 100 atoms. The approximate axial and 
radial dimensions of the atomic cloud are w z — 30 /im 
and w T = 2.8 /im, where w denotes the standard devi- 
ation of the density distribution. The signal photons, 
resonant with the 5s 2 S 1/2 (F = 2) 5p 2 P 3/2 (F = 3) 
transition in 87 Rb at 780.2 nm (see Fig. lb), propa- 
gate along z. The signal beam is focused to a 1/e 2 
radius of 1.2 ± 0.1 /im at the centre of the atomic en- 
semble. A counter-propagating control beam with wave- 
length 480 nm couples the signal transition to a highly 
excited Rydberg state with principal quantum number 
n = 60. The control beam is focused to a 1/e 2 radius 
of 17.9 ± 0.3 /im. The peak value of the control and 
signal beam Rabi frequencies are £lc/2ir = 6 MHz and 
f2s/27r = 1.2 MHz, respectively. For these parameters, 
the EIT width is A E it/2tt = 3 MHz and the blockade 
radius is R « 6 /im. 

In Fig. 2 wc illustrate the photon storage and retrieval 
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FIG. 1: Schematic of the experiment, a The atomic level scheme used to interface optical photons to the microwave field, b 
Optical photons are stored as Rydberg polaritons in a cold atomic ensemble. Subsequently, the internal states of the polaritons 
and their interactions are controlled using a microwave field. Finally, the modified optical field is read out and detected using 
time-resolved single photon counters arranged in a Hanbury Brown-Twiss configuration. 



process. The signal pulse is stored by reducing the inten- 
sity of the control field over a time of 100 ns. After the 
desired storage time, the control field is turned back on to 
read out the polariton field. This cycle is repeated every 
6 [as. The retrieved signal is typically around 200 ns long, 
which is determined by the control field switching time. 
The corresponding bandwidth of the storage process, A s , 
is 1.34±0.04 MHz (FWHM). There is some leakage of the 
polariton during the storage interval due to spontaneous 
emission and imperfect extinction of the control field. In 
addition, we note that the storage efficiency is far from 
optimised. However, efficiencies approaching 100% are in 
principle feasible by mode-matching to the time-reversed 
single-photon emission process [24|. 



Dipolc blockade limits the number of excitations that 
can be written into the sample leading to anti-bunching 
in the read out pulse as demonstrated in [ll[. To ob- 
serve this photon blockade effect, we perform a Hanbury 
Brown-Twiss coincidence measurement on the retrieved 
photons (see Fig. lb). The photon coincidences charac- 
terized by the second-order correlation </ 2 ) is plotted in 
Fig. 2c. There is a peak every 6 /is corresponding to the 
repetition rate of the experiment. In the absence of pho- 
ton interactions, the height of the peaks is expected to be 
unity indicating no bunching or anti-bunching (in prac- 
tice, variation in storage efficiency throughout an exper- 
imental run leads to a level 1.101 ± 0.003). In contrast 
the probability of coincidences within each pulse gives 
ff ( 2 )(0) = 0.78 ± 0.04. This partial suppression of g {2) (0) 
is consistent with a sample that is longer than the block- 
ade radius as shown in Fig. lb. 



Each photon is stored in the collective polariton state 
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and N is the number of 
. The phase factors are 



given by <pj = k ■ fj, where k is the effective wavevector 
of the spin wave, and r 3 - is the position of atom j. The 
phase of each term in the superposition ensures that the 
read-out emits a photon into the same spatial mode as 
the input. The lifetime of this phase- matched polariton 
is limited to of order microseconds by motional dephas- 
If atomic motion were reduced by additional 
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ing 

cooling, the decoherence time would be ultimately lim- 
ited by Rydberg lifetime which scales as n 3 . For 60sx/a 
the Rydberg lifetime is of order 100 /is. 

Coherent control of the stored photon is performed us- 
ing a resonant microwave field to couple the initial col- 
lective state Is) to a collective state 
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where p denotes an np Rydberg excitation. The states 
|s) and \p) form a two- level basis for collective encod- 
ing of the stored photon (24|. As the dipole moment for 
the ns —> np transition scales as n 2 , the figure-of- merit 
for single qubit rotations (Rabi frequency x dephasing 
time) scales as n 5 . For n = 60, of order 1000 qubit ro- 
tations within the decoherence time are possible. In this 
collective basis, both the Rabi oscillation frequency and 
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FIG. 2: Photon storage and retrieval: a The photon storage 
process begins at t = when the control field (blue) is turned 
off. At this time optical photons from the signal field are 
stored as Rydberg polaritons. The signal pulse (red) termi- 
nates at t — 0.6 fis. After a storage time of roughly 900 ns the 
control field is turned back on to read out the polariton field. 
The retrieved signal (filled circles) appears as a peak with a 
FWHM of 120 ±20 ns. The background signal without atoms 
(open circles) is shown for reference. A black band highlights 
the time window taken as the retrieved signal, b The detuning 
dependence of the retrieved signal indicates the bandwidth of 
the storage process. The retrieved signal is normalised by the 
total number of store/retrieve experimental cycles. The data 
are the mean ± standard error of ~ 2 x 10 J cycles per point. 
The line is a Lorentzian fit with FWHM of 1.34 ± 0.04 MHz 
(see main text), c The normalised second-order correlation 
function of the retrieved signal, binned over the entire 
retrieved pulse, as a function of the time delay r between the 
two detectors. The suppression at r = is a signature of 
dipole blockade during the polariton write process. 



the dephasing rate are independent of the atom number 
N (in contrast to the transition from the ground state 
1 0i • • • On) to the collective state |s), where the Rabi fre- 
quency scales as y/~N [IB])- This is important, as it allows 
us to observe the high contrast oscillations over many cy- 
cles even for a non-deterministically loaded sample. 

We study Rabi oscillations for n — 60 polari- 
tons coupled to a microwave field resonant with the 
6OS1/2 — >59p3/2 transition at 18.5 GHz. As the read out 
state is |s), the retrieved signal oscillates between a max- 
imum when the polariton state is \s) and a minimum 
when the polariton state is \p). In Fig. 3 we plot the 
retrieved photon signal as a function of the rotation an- 
gle, = fl^t. The microwave pulse duration is fixed so 
the microwave Rabi frequency f2 M increases from left to 



right in Fig. 3. Counter-intuitively, the Rabi oscillations 
revive for large 0. 

To understand these unusual dynamics, consider the 
pairwise dipole-dipole interaction between the cycling 
Rydberg polaritons (see also the Supplementary Infor- 
mation). The resonant dipole-dipole interaction energy 
between polariton modes is Vdd — d 2 /(4ireoR^), where 
R Q is the correlation length associated with the 60s block- 
ade process. For the range of Rabi frequencies shown in 
Fig. 3a we change between a regime where fi^ < Vdd 
on the left hand side and f2 M > Vdd on the right hand 
side. For f2 M < Vdd, the resonant dipole-dipole inter- 
action associated with the microwave transition domi- 
nates. In this case, the blockade sphere associated with 
the microwave transition is larger than the blockade 
sphere associated with the formation of the 60s polari- 
tons, R^ > R Q . As the resonant dipole-dipole interac- 
tion is an exchange process, this regime is dominated by 
excitations hoppin g, le ading to loss or dephasing of the 
polariton read-out [34[ . Consequently the retrieved pho- 
ton signal is suppressed and fits to an exponential decay. 
In this dephasing regime one may expect only a single ex- 
citation to survive and hence strong anti-bunching in the 
retrieved photon signal (34|. This effect was not observ- 
able in the current experiment, as the leakage apparent 
in Fig. 2a contributes a larger fraction of the signal in 
the dephasing regime and dominates the photon statis- 
tics. In contrast, for f2 M > Vdd strong driving forces 
the dipolcs to oscillate in phase, which suppresses the 
out-of-phase exchange interaction. In this case, where 
the microwave blockade radius is smaller than the opti- 
cal blockade, R^ < R Q , the exchange or hopping term is 
reduced to V d \/Q,^. Consequently the spin wave dephas- 
ing is reduced and the Rabi oscillations reappear. This 
recovery in the Rabi oscillations is a direct signature of 
the spatial correlations between Rydberg polaritons, and 
occurs when the microwave Rabi frequency is sufficient 
to overcome the polariton-polariton blockade. 

Significantly, A/"-particle correlations in the read out 
give rise to enhanced sensitivity to the rotation angle 0, 
which could be exploited in quantum metrology applica- 
tions [l6[. The retrieval probability, 



P 



is given by applying a Wigner rotation matrix to the 
collective Dicke state of N spins (see section I of Sup- 
plementary Information). This many-body character of 
the collective read-out is clearly visible in Fig. 3b. If 
we fit to [cos 2 (fi M t/2)] A ^ allowing M to float we obtain 
Af = 3.0 ± 0.2. This is consistent with the number of 
blockade spheres in our ensemble, given the 60s blockade 
radius and the geometry of the atomic cloud, see Fig. lb. 

The ability to control the quantum state of stored 
photons opens some interesting prospects for advances 
in quantum information and quantum simulation of 
strongly correlated systems. For example, the com- 
petition between resonant energy exchange (hopping) 



4 




FIG. 3: Controlling the interaction between Rydberg polaritons. a The retrieved signal, normalised to the case where no 
microwave coupling is applied, is plotted as a function of the microwave Rabi frequency, 57 M . The dynamics depend on the ratio 
between the Rabi coupling and the dipole-dipole interaction f2 M /Vdd - the condition £1^ — Vdd is indicated by the dashed line. 
For S7 M < Vdd (left hand side), resonant energy exchange between polaritons dominates over Rabi oscillations. This destroys the 
phase of the stored spin wave which is observed as a suppression of the retrieved photon signal. For fl M > Vdd, Rabi oscillations 
dominate and the exchange process is suppressed as the strong driving lifts the degeneracy between the dipole-dipole coupled 
states. The solid line is a phenomenological fit using the characteristic form for A/-particle Rabi oscillations coupled to a single 
optical read out mode, P = [cos 2 (n fi t/2)]' v . This function is combined with a tanh envelope, and an exponential decay at low 
microwave Rabi frequencies. From the fit we obtain M = 3.0 ± 0.2. Inset: Spin model of the dynamics. The dipole-dipole 
interaction (circles between atoms) favours excitation exchange between out-of-phase atomic spins (straight arrows) whereas 
strong microwave driving (circles around atoms) favours in-phase oscillations and suppresses the exchange process, b Higher 
resolution data of Rabi oscillations in the strong driving regime. The data are the mean ± one standard error of ~ 10 5 
store/retrieve cycles per point. The line is a similar fit to Fig. 3a, with M = 3.0 ± 0.2. 



and localisation is reminiscient of the Jaynes-Cummings- 
Hubbard model [13] (see section II of Supplementary In- 
formation). The ability to control the interactions be- 
tween polaritons using microwave fields provides a viable 
route towards fully deterministic photonic phase gates us- 
ing single photons , or to generate non-classical states 
of light from classical input fields [28[ . Finally, Rydberg 
polaritons provide a convenient interface between quan- 
tum systems that operate in the microwave and optical 
domains, such as circuit QED Q and atomic ensembles, 
respectively. Rydberg polaritons act as a source of quan- 
tum light, that can be coupled to on-chip [29| microwave 
resonators which in turn interface to solid state qubits 
@, H ' forming a complete architecture for transmitting, 
storing and processing quantum information. 



Methods Summary: Laser cooled 8 Rb atoms are 
loaded into an optical dipole trap at A = 910 nm. The 
trapping light, which co-propagates with the signal beam, 
is focussed to a 1 /e 2 radius of approximately 5 /xm. The 
trap depth is roughly 0.6 mK and the temperature of the 
trapped atomic cloud is typically w 100 /iK. After rc- 
pumping the atoms into the 5si/2(-F = 2, ray = 2) state, 
the dipole trap power is switched on-and-off with a pe- 
riod of 6 fis and a 50% duty cycle. Experiments are per- 
formed during the periods when the trap is off to avoid 
differential light shifts, with one store/retrieve experi- 
ment per trap-off period. The trap modulation typically 
lasts 20 ms, after which the storage efficiency begins to 



decrease. The trap is then reloaded and the sequence is 
repeated. 

When the dipole trap has been switched off, a signal 
pulse of 1.1 ^s duration is stored by reducing the intensity 
of the control field. A microwave pulse of duration 150 — 
300 ns is applied, coupling the 60s!/2 and 59p 3 / 2 Rydberg 
states. The dipole moment for this transition is d = 
\/2/9 x 3468 eao, where the radial matrix element was 
calculated using the Numerov method. The microwaves 
are emitted from a stub antenna, orientated such that the 
microwaves drive 7r-transitions. Following the microwave 
pulse, the control field is turned back on to read out 
the polariton field. Rabi oscillations are mapped out by 
varying the power of the microwave pulse, keeping the 
pulse length fixed. The total storage time of the polariton 
is approximately 900 ns. 

The retrieved photon signal is windowed over the re- 
gion highlighted in Fig. 2a. When calculating g^ 2 \ only 
photons within this time window are correlated. For the 
Rabi oscillation measurements in Fig. 3, the total num- 
ber of photon counts within the window is extracted for 
each shot. There are 30 shots for each data point, with 
3334 individual experiments in each shot. 
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SUPPLEMENTARY INFORMATION 



A. Collective Rabi oscillations of J\f indistinguishable Rydberg polaritons 

Photon storage in the regime of dipole blockade forms the polariton superposition state 

N 



TV . 

where |s J ) = IO1O2 . . . sj . . . Oat) is a state with the single excitation at atom j within an ensemble of N atoms. The 
spatial dependence of the phase factors <f>j = k ■ fj, where k = k s + k c , is determined by the sum of the wavevectors of 
the excitation lasers and the position of each atom fj. The phase factors determine the directionality of the optical 
read out of the polariton. 

A resonant microwave field couples the singly-excited, symmetric polariton state |s) to a similar state 



1 N 

\p) = 77jjY,^ 3 \P°) ■ ( 2 ) 



i=i 



Under microwave coupling each Rydberg polariton (neglecting polariton-polariton interactions) becomes an effective 
two-level system or spin-^ quasi-particle. As the wavelength of the microwaves is much larger than the sample size, 
the global phase-structure of the polariton is preserved. Any operator acting on this system can be written in terms 
of the Pauli spin matrices a x , a v , and u z . For Af polaritons we define collective spin operators 



(3) 



and similarly for y and z. The Dicke states labelled \J,M) are the eigenstates of J 2 = J 2 + J 2 + J 2 and J z . For 
coherent driving, the evolution of the A/"-particlc Dicke states is given by the Wigner rotation matrix T>. The reduced 
rotation matrix between states with the same J but projections M' and M is [301 ] 



& (oeo) - I (J-m)\(j + m')\ cos2J+ m-m> f®\ sin M'-M f®\ p U) 

V M , M {0, 0, 0) - (M , _ M)! W (J + M)|(J _ M/)! cos ^ 2 J sin [ 2 J ^ (4) 



where F = 2F\(M' — J, —AI— J; M' — M + l; — tan 2 8/2) is a Gauss hypergeometric function, and the angle = fl^t 
is the rotation induced by the microwave field. The initial spin state is | J = AT/2, M — —Af/2). To retrieve the stored 
optical pulse in the same mode we need to conserve the phase-structure, therefore we project onto the identical state 
|J = Af/2, M' = — Af/2). In this particular case, where M = —J, the hypergeometric function takes the value 1. The 
probability of photon retrieval is proportional to the matrix element squared giving 



P = 



cos 2 



2 



(5) 



This result gives an accurate description of the collective read out of the rotation of the A/"-polariton state. Note 
that the frequency of the oscillations corresponds to the single-atom Rabi frequency , and does not depend on the 
number of polaritons Af, or the number of atoms in each polariton N (assuming that TV is sufficiently large to ensure 
directional emission). This feature is important, since it allows us to observe the oscillations in a non-detcrministically 
loaded sample. The power of the cosine depends only on the mean number of excitations Af arising during to the 
photon storage protocol. Larger Af gives enhanced sensitivity to the rotation angle around integer multiples of 2tt. 

Finally, it is important to point out that this picture only applies if the number of atoms, N, in each polariton, and 
their spatial extent are sufficiently large to support directional emission during read out. The condition on the size 
is that the blockade radius should be many times the wavelength of the optical transition R a 3> A. Both the size and 
number conditions break down at low principal quantum number, n, and the dynamics revert to single particle Rabi 
oscillations. 
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B. Spin chain of strongly interacting polaritons 

If we treat the Rydberg polaritons as localised spin-i quasi-particles with a separation R Q , then the interaction 
Hamiltonian for zero detuning and only nearest-neighbour interactions is 

TV TV 

H = a a < + + Kid °t°~i + f E( a * + a l) > ( 6 ) 

t=l (ij) i=l 

where the first term is the usual Jaynes-Cummings coupling between a two level system and a quantised microwave 
field, the second term with Vdd = d 2 /4ne Rl describes the dipole-dipole interaction between sites i and j, and the 
third term describes coherent driving. We note in passing the analogy to the Jaynes-Cummings Hubbard Model 
(JCHM) which describes photon tunnelling between arrays of cavities with one two-level system in each cavity [3lj |. 
The main difference is that whereas in the JCHM real photons hop between cavities, for Rydberg polaritons there is 
an exchange of virtual photons, i.e., the photon hopping term a\a,j in the JCHM is replaced by an excitation hopping 
term erf <j~ . 

The current experiment differs from the situation described by © in a number of ways. First, the absence of a 
microwave cavity to enhance the single photon coupling g means that we can treat the microwave field classically. 
Second, the field drives a 7r-polarizcd sj/2 — > P3/2 transition but the dipole-dipole interaction can also couple to other 
angular momentum states. Neglecting the fine and hyperfine structure, the interaction Hamiltonian has the form 

n N 

h = + Vdd I>i>7 + ^4 - W) . ( 7 ) 

i=X (ij) 

where is the Rabi frequency on the microwave transition and [i i ' is the dipole operator for Am = 0, ±1 transitions 
on the i-th polariton. The eigenvalues for two localised polaritons as a function of their separation are shown in Fig. 
SI. The energy scales can be parameterized in terms of the two radii, R D — (Cs/^Aeit) 1 ^ 6 and i? M = (C^/ftf?^) 1 / 3 , 
characterising the optical and microwave blockade, respectively. For 60s60s Cq = —140 GHz fim 6 and for 60s59p 3 / 2 
C3 = —14.3 GHz ^m 3 [H, HI]. For weak driving, Fig. Sla, there is a strong state mixing at R which suppresses 
the retrieved photon signal whereas for strong driving, Fig. Sib, the dressed states are only weakly perturbed by 
dipole-dipole interactions. 




Rlfjm Rljjm 

FIG. 4: Eigenvalues of the polariton Hamiltonian (7) as a function of the polariton spacing. The blockade radii for optical, 
Ro — 6 /J,m, and microwave, Rp,, excitation are indicated, a For weak driving — 20 MHz, the dipole-dipole interaction 

dominates, _R M > R , and there is a strong mixing of states with different angular momentum, m. b In contrast for strong 
driving fi M /27r = 200 MHz, _R M < R , the splittings between the coupled dressed states are only weakly perturbed by the 
dipole-dipole interaction (of order V^ d /Q^) and there are no level crossings for R > R . 

A final consideration is that the polaritons are delocaliscd. Consequently, the dipole-dipole exchange terms are 
different for each atom pair, /x, v, and we should replace Vdd with a sum over individual atoms in each polariton 
V^i ■ I n tne re gi me : > Ro, this gives rise to loss and dephasing of the spin wave [U as illustrated in Fig. 3a. 
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